New mechanism of membrane fusion 
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We have carried out Monte Carlo simulation of the fusion of bilayers of single chain amphiphiles 
which show phase behavior similar to that of biological lipids. The fusion mechanism we observe is 
very different from the "stalk" hypothesis. Stalks do form on the first stage of fusion, but they do 
not grow radially to form a hemifused state. Instead, stalk formation destabilizes the membranes 
and results in hole formation in the vicinity of the stalks. When holes in each bilayer nucleate 
spontaneously next to the same stalk, an incomplete fusion pore is formed. The fusion process is 
completed by propagation of the initial connection, the stalk, along the edges of the aligned holes. 
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I. INTRODUCTION 

Fusion of membranes is involved in basic biological 
processes but its mechanism remains poorly understood. 
That different proteins trigger different fusion events ob- 
scures the possibility of a stage common to them, one 
which depends only on the properties of membrane bi- 
layers themselves. The existence of such a universality is 
made plausible by the realization that, in all fusion pro- 
cesses, membranes must merge, and the properties of the 
membranes will most likely determine the nature of this 
stage, the fusion intermediate.!!! Time and length gcales 
of fusion evejats are of the order of microsecondsH and 
nanometers,™ so that direct measurements of intermedi- 
ate structure have not been possible. Limited theoretical 
treatment of this problem, based on membrane elasticitsz. 
theory, has focused on the so-called stalk mechanism.EllD 
In this scenario, the leaves of the two cis membranes, 
closest to one another, fuse forming a stalk. The stalk 
expands radially and thins so that the trans leaves make 
contact and form a single bilayer. This stage is denoted 
hemifusion. Hole formation in this bilayer completes the 
fusion pore. Due to the lack of direct experimental con- 
firmation, this mechanism, although plausible, remains 
hypothetical. 

To obtain a direct view of the fusion mechanism we 
have undertaken the first Monte Carlo simulation of 
membrane fusion, one in which membranes are treated 
on the molecular level. In contrast to the stalk hypothe- 
sis, we find a very different fusion mechanism. Although 
stalk formation does take place in the initial stage of 
fusion, it does not result directly in the fusion pore. In- 
stead, the stalk destabilizes the contacting bilayer mem- 
branes and promotes formation of transient holes next to 
it. When two holes are nucleated in apposing bilayers in 
the vicinity of the stalk, the latter grows along the hole 
edges, sealing their periphery like a zipper, and eventu- 
ally closes thereby forming the fusion pore. 



II. THE MODEL 

Our model describes membranes formed by single- 
chain amphiphiles, like block copolymers which exhibit 
the sarnfi phases as do biological lipids, and also form 
vesicles.! It has the advantage that it has been well stud- 
ied, permits detailed analysis of molecular configurations, 
and is well suited to processes occurring on the small 
time and length scales characteristic of fusion. The arru 
phiphiles are treated using the bond fluctuation modelt! 
in which each molecular segment occupies a cube of a 
three-dimensional lattice. The eight lattice sites defin- 
ing the cube cannot be occupied by another segment 
centered on neighboring sites. Segments along an am- 
phiphile are connected by one of 108 bond vectors of 
lengths 2, VS, V5, 3 or \/T0, measured in units of the 
lattice spacing oq. Mapping this model onto lipids in 
solutionJiJ we find the lattice spacing to correspond to 
approximately lA. The large number of bond vectors 
and the extended segment shape allow a rather faithful 
approximation of continuous space, while retaining the 
computational advantages of lattice models. The am- 
phiphilic molecules consist of = 32 segments, of which 
10 are hydrophilic and 22 are hydrophobic. This partic- 
ular choice of the ratio of hydrophilic and hydrophobic 
segments results in the diblock system being close to co- 
existence of the lamellar and inverted hexagonal phases. 
The solvent is represented by a homopolymer, chains con- 
sisting of 32 hydrophilic segments. Like segments at- 
tract each other and unlike segments repel each other 
via a square well potential which comprises the nearest 
54 lattice sites. Each contact changes the energy by an 
amount e = QAllksT. The particular choice of the inter- 
action parameter e guarantees that the interfacial width 
between hydrophilic and hydrophobic segments is not too 
small to be comparable to the lattice spacing, and at the 
same time results in well-defined bilayers. 

The simulation cell is L x L in the x,y directions and 
of length D in the z direction, with L — 156ao and 
D = 96ao. Periodic boundary conditions are utilized 
in all three directions. The monomer density of the sys- 
tem is p = l/(16ao), corresponding to 146,016 segments 
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within the volume, or 2376 amphiphiles and 2187 ho- 
mopolymers. To encourage the fusion of the bilayers that 
we study, we prepare them under tension by providing 
only enough amphiphiles to form bilayers of thickness 
25ao, thinner than their equilibrium thickness of 25.2ao, 
which we determined independently. Given that the sys- 
tem is nearly incompressible, this corresponds to a frac- 
tional increase of area of less tha£L|£L8%, one easily sus- 
tained by biological membranes. Ilirlij Two such bilayers 
are created parallel to the x — y plane, stacked one upon 
the other, thereby mimicking the dehydration which per- 
mits close bilayer contact, a circumstance known to pro- 
mote fusion.Ej Sixty-four independent starting configura- 
tions were prepared. Monte Carlo simulations were per- 
formed in the canonical ensemble. The conformations are 
updated by local monomer displacements and slithering- 
snake like movements. The different moves are applied 
with a ratio 1 : 3. The latter moves do not mimic the re- 
alistic dynamics of lipid molecules and we cannot identify 
the number of Monte Carlo steps with time. However, 
the density of hydrophilic and hydrophobic segments is 
conserved and the molecules diffuse. Consequently, we 
expect the time sequence on length scales larger than a 
single molecule to resemble qualitatively those of a sim- 
ulation with more realistic dynamics. 

To make sure that an isolated membrane is stable, we 
have simulated a single bilayer under the above condi- 
tions. The system was extremely stable with respect to 
hole formation, apparently due to a very high energy cost 
to form the hole's edge, and exhibited the usual capillary 
wave fluctuations. In contrast, the system with two ap- 
posing membranes resulted in a range of structural trans- 
formations that eventually led to fusion pore formation. 

III. RESULTS 
A. Qualitative observations 

Initially the apposed bilayers formed a number of lo- 
cal connections, the so-called stalks. Formation of these 
stalks was promoted by the system being near bulk coex- 
istence between lamellar and inverted hexagonal phases. 
In the stalk mechanism, these stalks grow radially and 
eventually form transmembrane contacts. In contrast, 
we observed anisotropic growth of the stalks, which re- 
sulted in the formation of structures whose length was a 
few times that of their width. These may be precursors 
of line defects observed during the lamellar to inverted- 
hexagonal transition.!!! 

Formation of the elongated stalks led to local destabi- 
lization of the bilayers. Holes formed next to these de- 
fects. We stress again that single, isolated, membranes 
were very stable. Thus hole formation in the apposed bi- 
layers is apparently due to inter-bilayer interactions ex- 
pressed in the stalk formation. It is easily seen that a 
hole formed next to a stalk has a smaller edge energy, in 
comparison to an isolated hole due to the reduction in 
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FIG. 1. The fusion intermediate, (a) Schematic sketch of 
the intermediate. Only hydrophobic portions are shown, (b) 
Top view of a fusion pore through two membranes. Only 
hydrophobic segments are shown, (c) Top view of the layer 
between the two membranes. Hydrophilic segments are dark 
grey, hydrophobic segments are light grey, (d) Side view of 
the snapshot. 

curvature energy. This reduction is roughly proportional 
to the length of the edge adjacent to the stalk. 

Hole formation was a transient process, i.e. numerous 
holes were observed to open and close reversibly. Only 
when two holes, one in each bilayer, nucleated next to, 
and on the same side of, a given stalk did the final stage of 
the fusion process commence. The stalk then propagated 
along the edges of the two holes bringing them together, 
like a zipper, to complete the fusion pore. 

We show what one such partially formed pore looks 
like in FIG. |l|. Panel (b) presents the top view of the 
two membranes. Only tail segments are shown. The 
pore through the two bilayers is clearly visible. Panel 
(c) presents the top view of the layer between the two 
membranes. Hydrophobic segments are light grey, while 
hydrophilic segments are dark grey. Connectivity is es- 
tablished along a portion of the pore's rim by the stalk, 
while the two hydrophobic sheets are still unconnected 
elsewhere. This is confirmed by a side view through the 
pore (d). A schematic of this intermediate is shown in 
panel (a). 

B. Quantitative analysis 

After 62,500 attempted local monomer displacements 
per segment and 187,500 slithering-snake movements per 
molecule, we find that the initially sharp interface of the 
bilayers have widened, and profiles across the membranes 
have adopted locally their stationary form. Profiles at a 
later stage are broadened due to fiuctuations of the local 
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membrane position, and to the formation of holes in the 
membranes. In order to reduce the effect of membrane 
fluctuations on our observation of the local profiles, we 
define at each x and y a local midpoint, 2:middio(a;, ?/), of 
the two-bilayer complex as follows .llJ For each x, y and z, 
we determine the total number of amphiphile segments, 
na{x,y, z), in a volume centered at these coordinates, of 
width i? = 28 oo in the lateral directions, and 50 aq verti- 
cally, the thickness of two bilayers. The value of the coor- 
dinate z which maximizes na{x, y, z) defines Zmiddic(2;, y). 
We have chosen a rather large lateral length scale, B, so 
that a small hole in a single bilayer does not significantly 
alter our estimate for the midpoint of the two bilayers. 

To analyze structural changes in the system, we looked 
not only at the fusion pores traversing both bilayers, but 
also at the holes in each bilayer which we found to occur 
about an order of magnitude more frequently than pores. 
We did this by calculating the local density of hydropho- 
bic tail segments coarse grained over 4ao x 4ao square 
columns centered at x and y. Holes and pores should 
have a very low, or vanishing, density of tails compared 
to the areal density of two normal bilayers. If the local 
value of the coarse-grained density of tail segments within 
a column extending through the whole z-range falls be- 
low 25% of the normal areal density of the two bilayers, 
we define this column as belonging to a pore traversing 
them. To locate a hole in only one of the bilayers, the 
integration over the z coordinate, rather than extending 
over the whole simulation box width, D, is limited to ci- 
ther the upper or lower half of it starting at the point 
between the two bilayers. If the local coarse-grained tail 
density falls below 25% of the normal areal density of 
a single bilayer, the column belongs to a hole. Pores 
and holes are defined as the aggregate of nearest neigh- 
bor 4ao X 4ao plaquettes with low hydrophobic segment 
density. 

The probability of hole sizes in the two apposed bi- 
layers and the isolated bilayer are presented in FIG. ||. 
Not only is the absolute frequency of hole formation very 
much less in the isolated membrane than in the two ap- 
posed bilayers, but its decay with increasing hole size is 
also much more rapid. This indicates that the edge free 
energy of holes in the isolated bilayer is larger and cor- 
roborates the observation that stalk formation promotes 
hole formation when two bilayers are in contact. 

To determine unambiguously whether fusion pore for- 
mation involves a stage in which the trans leaves make 
contact, one unavoidable in the standard stalk mecha- 
nism, we perform a statistical analysis of the fusion in- 
termediates observed in our simulations. We examine in 
each simulation run the first pore which exceeds a cer- 
tain radius, Rmm = 5.5ao, one sufficiently large that the 
pore must be open; that is, there is a channel connecting 
upper and lower parts of the system accessible to solvent 
molecules. 

Our strategy is to look at the region along the rim, 
or edge, of the pore in the pore's mid plane. If fusion 
proceeds via hemifusion this region has to be essentially 
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FIG. 2. Normalized probability distribution of the hole 
sizes in a bilayer. Circles refer to the system of two apposed 
bilayers, while squares denote the results for a single bilayer. 
Data for the two apposed bilayers were collected from 64 sys- 
tems and 750000 up to 956250 slithering-snake steps. Data 
for the single bilayer were sampled over 7 configurations but 
much longer times. 



completely hydrophobic. In contrast, if fusion proceeds 
through the formation of two transient holes in apposing 
bilayers, then we may see pores in an intermediate stage 
in which the holes are aligned, but the rim of the pore has 
not fully formed. In such a case, the region along the rim 
of the pore would be at least partially hydrophilic. We 
measure the density of segments in a series of columns of 
height 5ao and base 4ao x 4ao centered along the rim of 
the pore, which is defined as the points in the x — y-mid 
plane with a distance 14ao ± 2ao from the center of the 
pore. This distance corresponds to a sum of radius of the 
pore, 5.5ao, and half of the hydrophobic thickness of the 
bilayer, 8.5ao. 

To quantify the connectivity of the pore's rim, we ex- 
tract the distribution of the density difference A(f> be- 
tween hydrophilic and hydrophobic segments in this re- 
gion. On average there are = 5 monomers in a column 
whose volume is 4ao x 4ao x 5ao. If the connectivity of 
the hydrophobic regions of the two bilayers were com- 
plete, as in the stalk mechanism, then we expect the 
rim region to be hydrophobic, so that the distribution 
of A(j) from all these columns would exhibit a single peak 
around A0/0o = ^1- This is not at all what we see. In- 
stead, only 12 pores have this single peaked distribution, 
whereas the rest are either bimodal (37 configurations) 
with a second peak centered about A(j)/(f)o = -1-1 or have 
a peak at A0/0o = —1 with a long "tail" (15 configu- 
rations) toward positive A(f). From this we can conclude 
that in the majority of cases there is a considerable por- 
tion of the rim which is still hydrophilic; that is, the rims 
of the two holes are only partially connected. The aver- 
age of all distributions which we obtain are shown by the 
solid circles in FIG. ||. The distribution of the values of 
A(f) obtained from the series of columns along the rim of 
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FIG. 3. Probability distribution of the difference between 
hydrophilic and hydrophobic segment density in the rim. Pos- 
itive and negative values correspond to hydrophilic and hy- 
drophobic dominated regions respectively. Circles show the 
average over all 64 independent configurations, while the sohd 
line depicts the result for the snapshot in FIG. ^ 



the particular pore in FIG. |^ are shown as a solid line. 



IV. DISCUSSION 

The fusion mechanism W|e-.propose is similar to that 
considered for electrofusionjld with the addition that we 
specify the form of the fusion intermediate. It is in ac- 
cord with the recent observation in molecular dynam- 
ics simulations that lipid tails take on a much greater 
range of configurations, including those sampling the hy- 
drophilic-TCgion, when membranes are brought into close 
contact.EEl Our mechanism is similar in spirit to the initial 
stages of the mechanism of hemagglutinin-mediated virjal 
fusion recently proposed by Bonnafous and StegmannEZI 
who suggest that before fusion can occur, a hole must 
first form in the target membrane. However they go on 
to suggest that this hole leads to fusion of the cis lay- 
ers producing a hemifusion diaphragm in which a hole is 
nucleated. In contrast, hemifusion plays no role in the 
scenario we propose. We do observe strong mixing of am- 
phiphiles in the cis leaves before pore formation. From 
such mixing, it is-often inferred that hemifusion precedes 
complete fusion.Ea Instead, we find that this exchange 
results simply due to formation of the stalks. 

The destabilization of the membranes by the stalks 
that leads to formation of the transient holes in both bi- 
layers, observed-hy us, is consistent with the work of Cevc 
and Richardsenfij who emphasize the fact that membrane 
fusion is strongly promoted by defects in the bilayer struc- 
ture. It seems clear that the additional destabilization of 
bilayers by the defects embedded in the membranes prior 
to fusion facilitates formation of the transient holes nec- 
essary for fusion. 



The destabilizing role of the stalk intermediates, ob- 
served in our simulations, is missing from the current 
phenomenological models of membrane fusion. Although 
we do not observe the standard stalk mechanism, our sim- 
ulations do not rule it out under different thermodynamic 
conditions. Extensions of phenomenological approaches 
that include the new observed intermediates should clar- 
ify this point. We expect the fusion intermediate we have 
seen in our simulation to be similar to that in biological 
membranes. Its elucidation should facilitate the control 
and modification of the fusion process itself. 

Note: While preparing this manuscript, we learned 
that Noguchi & Takasu,Ej studying a very different sys- 
tem consisting of rigid amphiphilic molecules of three 
atoms and containing no solvent, observed fusion behav- 
ior over a limited temperature range which is similar to 
that we observed in our system. This may indicate a 
universality to the mechanism we have described. 
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